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Antiferromagnetic Hamiltonians with short-range, non-frustrating interactions are well-known to 
(^-^ , exhibit long range magnetic order in dimensions, d > 2 but exhibit only quasi long range order, with 

f*"*) ■ power law decay of correlations, in d = 1 (for half- integer spin) . On the other hand, non-frustrating 

' long range interactions can induce long range order in d — 1. We study Hamiltonians in which the 

f\| i long range interactions have an adjustable amplitude A, as well as an adjustable power-law l/l^l™, 

using a combination of quantum Monte Carlo and analytic methods: spin-wave, large- iV non-linear 
a model, and renormalization group methods. We map out the phase diagram in the X-a plane 
and study the nature of the critical line separating the phases with long range and quasi long range 
order. We find that this corresponds to a novel line of critical points with continuously varying 
critical exponents and a dynamical exponent, z < 1. 
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I. INTRODUCTION 

The ground-state (GS) of the nearest neighbor antiferromagnetic (AF) Heisenberg model on a bipartite lattice: 

h = £3-4 (i) 

(id) 



is generally expected to have long range order (LRO): 

a 

g; (So • S r ) - ±m 2 AF , (r->oo), (2) 

for any spin magnitude, S and any dimension d > 2 On the other hand, in dimension d — 1, the behavior depends 
£SJ ' on whether S is integer or half- integer |2J. In the half-integer case the spin-spin correlation function 

o 

G\ ' (So ■ S r ) oc v ^' (3) 

r 

m . 

is expected characteristic of a guasi-long range order (QLRO). (In the integer spin, Haldane gap case, correlations 
decay exponentially.) This behavior in d = 1 for half-integer S is believed to be universal, not depending on the 
magnitude of S nor on the details of the Hamiltonian as long as it is short range and not too frustrating. Long range 
interactions (i.e. power-law decaying with the relative distance between interacting moments) can be introduced 
in spin models either for some experimental reasons like dipolar or RKKY interactions or simply because of some 
theoretical relevance. For instance a famous example is the Haldane-Shastry model with AF frustrating 1/r 2 
I ■ interaction which exhibits an exact RVB GS. Another theoretical interest comes from the possibility to interpolate 
between discrete dimensions by tuning continuously the exponent that governs the decay of the interaction with the 
distance. Indeed, the possibility for true LRO to occur in d = 1 with long range interactions has motivated many 
studies during the last decades 0, IS 0> @> II \Ei EH EH EE an d is the subject of the present paper. While a 
long standing debate about the critical behavior of the Ising model in d dimensions with long range ferromagnetic 
interaction decaying like r~ d ~ a has been quite intense during the last thirty years 0, ITTl Il2j the N-vector model 
has also been a subject of interest for many authors 0BISE9- Concerning the Heisenberg model with long range 
interaction ~ r~ a , the seminal paper of Mermin and Wagner [5| proving the absence of LRO at finite temperature T 
in d < 2 for a > d + 2 has been recently reconsidered by Bruno |13( who gave stronger conditions for the absence of 
spontaneous magnetic order at T > in d < 2. For instance, he proved that the AF non frustrating one dimensional 
model 



St-s^-x^i-iy^+i 

3=2 3 



(4) 



does not have Neel order for any temperature T > if a > 2. Actually, much less is known about the T = case, 
except the work of Parreira et al. |l5j where the authors signaled the existence of the bound a = 3 over which T = 
LRO is ruled out ^(|. A particular case, the model with A = 1 was recently analyzed at both T — and T > 



in Ref. [Tj], using the lowest order spin- wave (SW) approximation, expected to be valid for large enough S and small 
enough a. There it was shown that the SW dispersion relation takes the sublinear form, at low k: 



0) oc Ifcl^- 1 )/ 2 , 
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for a < 3. Consequently the quantum 1/5 reduction of the order parameter: 
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is finite for any a < 3. By requiring that Am, < S, a consistency condition on the SW approximation, it is concluded 
that LRO occurs for any S at sufficiently small a. (However, such an estimate is presumably only reliable for S 3> 1.) 
After correcting a numerical error in Ref. 0], the SW prediction for the S = 1/2, A = 1 case is existence of Neel 
order at T = for a < a s c w = 2.46. 

In this work, we extend the results of Yusuf et al. in several ways, focusing on the zero temperature behavior of 
the non frustrating spin 1/2 Hamiltonian (0} with long range interaction of adjustable strength A and exponent a. 
In Sec. II, we consider the relevance of the long range term as a perturbation to the nearest neighbor interaction, 
using a simple heuristic argument of mean-field type as well as the power-counting of the scaling dimension of the 
perturbation. For A < 1, we find that the long range perturbation is marginal if a = 2 and relevant (irrelevant) 
for a < 2 (a > 2). We then investigate the a- and A-dependence of the critical behavior using various techniques. 
We begin, in Sec. Ill, with semi-classical calculations: the SW expansion and a large-N approximation based on the 
non-linear a model. Both approximations give qualitatively similar phase boundaries, and sublinear dispersion like 
in Eq. (JSJ in the ordered phase. Some of the critical exponents can also be estimated within these approximations. 
However, the results obtained in the SW or large-N approximations are not quantitatively correct. We therefore use 
large scale numerical simulations to investigate more precisely the phase diagram of this model in Sees. IV and V. 
We study systems of up to L = 4000 sites usin g qu antum Monte Carlo (QMC) methods, based on a stochastic series 
expansion (SSE) of the partition function [rH Il8|. We verify that for S — 1/2, there are indeed stable phases with 
both QLRO given by Eq. 10 and with true Neel LRO [Eq. JSJ]. We accurately determine the phase boundary, as 
well as some of the critical exponents which are found to vary continuously along the critical line. In Sec. VI, we also 
apply analytic renormalization group (RG) methods to investigate the case A <C 1. Sec. VII contains conclusions. In 
two appendices we gives further details on the spin- wave theory and large-N calculations. 

II. RELEVANCE OF THE PERTURBATION: MEAN FIELD AND SCALING ARGUMENTS 

Let us consider a short range spin 1/2 chain with an additional long range perturbation of the form 



Following an argument given by Cardy [19j for the relevance of a long range perturbation, we can in first approximation 
look at the mean field correction to the free energy coming from this long range term (jTJ): 




(7) 



with 



J(r,r') 




(8) 
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where (...) is evaluated in the unperturbed system where we know the behavior of the correlation function 





(10) 
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In a finite system of length L, the change in the free energy per site Sf thus scales like 



3 



The integral above will give a constant term and a size dependent term 

5f(L) ~ i^-a-z--q ^ L~ a , (12) 

where we have used the fact that z = r\ = 1 in the short range QLRO regime of the spin 1/2 chain. Then, we can 
compare this with the usual finite size corrections to the free energy of the conformally invariant short range S — 1/2 
chain which are known to scale like L~ 2 to lowest order [2(j. This tells us that (to first order perturbation) if a < 2 
the long range perturbation creates a correction which dominates the L~ 2 correction of the unperturbed fixed point 
and is probably a relevant perturbation for the short range model. 

Another way of deriving this result is to compute the scaling dimension of the perturbation, based on the usual 
continuum formulation of the short range model in which uniform and staggered magnetization density operators, 
( Jl + Jr.) and n are introduced: 

(J L + J R ) + (-l) x n(x). (13) 

Only slowly varying Fourier modes of the fields Jl/r{x) an d n(x) are present in the low energy effective Hamiltonian. 

Jl/r ar e the conserved left/right-moving spin densities. Ignoring a marginally irrelevant interaction, the staggered 
magnetization field, n, has the Green's function: 

K(z)n"(0)) = f^, (14) 

with z = t + ix. The long range perturbation adds to the low energy, continuum limit of the imaginary time action, 
a term of the form: 

/di~ dx dy 
j——^n(T,x)-n(T,y). (15) 

Utilizing the fact that, from Eq. (|14|) . n has a scaling dimension of 1/2, a simple power counting tells us that the 
perturbation is irrelevant for a > 2, relevant for a < 2, and marginal for a = 2. Also note that A > corresponds to 
non- frustrating interactions which favor the Neel state with (n z ) ^ 0. 

III. SPIN- WAVE EXPANSION AND LARGE-N APPROXIMATION 
A. Spin-wave expansion 

This calculation simply generalizes that of Yusuf et al. in Ref. 0], to A ^ 1. We summarize here the main steps. 
Some further results are given in Appendix A. On the LRO side of the transition, we use the Holstein-Primakoff 
approximation [2lJ: 

S* = S — a|ai, 5+ w V^Sdi, S7 ~ V2Sa\, 

for i odd and 

S* - b]bj - S, S+ w V2Sb], Sr « y/2Sbj, 

for j even, and retain only the quadratic terms in the Hamiltonian ijljl. After a Fourier transform over the reduced 
Brillouin zone k E (—■£-,■£-), we find: 

n S w~Sj2[(l-f(k)) (a{a k + b\b k ) + g(k) (a\bl k + 6_ fc a fc )] + ... (16) 
k 

where, for an infinite chain 221: 
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l(k) = 2cos(fca) + 2A^ 



cos[fc(2n — l)a] 
(2n - l) a 



This quadratic Hamiltonian can be diagonalized with a Bogoliubov transformation to: 

T-Lsw ~ S^ujk [xt,iXk,i + xl, 2 Xk,2 / 

k 

with a SW spectrum 



"k = Vt7" /(*)]+ [ff(*)] S 
as discussed above. At T = 0, the correction to the staggered magnetization at any site is 

~7-f(k) 



(17) 



(18) 



Am q = (ajai) = (btbj) = — 



dk 



The consistency condition Am, < S then allows us to find the SW approximation for the value of af below which 
long range Neel order is established. As already stated, for 5=1/2 and A = 1, we find af" = 2.46. A plot of a^ w vs. 
A, for S = 1/2, is shown in Fig. ^ 




FIG. 1: Spin-wave approximation prediction for the value a™ below which long range Neel order is expected at T — 0, as a 
function of A and for S = 1/2. The critical curve saturates at a ~ 2.9032 when A — * oo. 



This phase boundary given by the lowest order SW approximation turns out not to be quantitatively correct, as 
we are going to see with the QMC results presented in section Hvl In particular it happens that SW predictions miss 
the fact that the critical line goes to a — 2 when A — > 0. Nevertheless SW predictions are, to some extent, shared by 
large- N calculations as we are going to see below. 



B. Large- N approximation 

The details of this calculation are given in Appendix [B] Here we present the main steps and discuss the results 
which come from this approximation. We generalize the Neel order parameter field, n(r, x) appearing in Eq. I|13|) . to 
an iV-component field and take the limit of large N . In this approximation 2 phases occur in the A -a plane. The 
critical line terminates at a = 1, as in spin- wave approximation. These two phases are a phase with Neel order and 
a disordered phase with a finite correlation length. (The unusual quasi long range ordered phase is special to the 
case N = 3 and is not captured by the large- N approximation.) Along the critical line separating these two phases 
the mean field result, rj — 3 — a is obtained. The dynamical exponent takes the value z = (a — l)/2, corresponding 
to the dispersion relation lu oc 

|fc|(a-i)/2 a i so obtained in spin-wave theory. The correlation length diverges with an 

exponent, v defined by: 



£cx|A c -Ar, 



(19) 



5 



with 



v = 1/(q - 1), (1 < a < 5/3) 
v = 2/(3 -a), (5/3 < a < 3). 



(20) 



IV. QUANTUM MONTE CARLO RESULTS I: FINITE SIZE EFFECTS 

In this section, we present results obtained using the QMC SSE method based on directed loop updates [l7j . 
This algorithm, used here to investigate the model Q, has been proposed recently by Sandvik I8| to study spin 
Hamiltonians with non-frustrating long range interactions. 



Finite size corrections 



We first focus on the A = 1 case, studied by SW in Ref. > governed by the following Hamiltonian 



t>J>l 



(-l)J ? a 



(21) 



In order to detect a Neel instability at the thermodynamic limit, we compute the staggered structure factor, normalized 
per site, on finite length spin S = 1/2 chains, defined by 



L 2 



(22) 



We have performed SSE simulations for different system sizes, up to L max — 4096, at temperatures /3 _1 = 1/2L low 
enough to get the GS properties. Results for S^L) are shown vs 1/L in the left panel of Fig. |3 for different values 
of the power-law exponent a. The staggered structure factor displays two types of behavior: for small values of a 
it saturates to a finite non-zero number whereas for large enough a, S W (L) vanishes when L — > oo. Then, in order 
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FIG. 2: Left panel: T — QMC results for the staggered structure factor per site S n (L) computed in the GS for the Hamiltonian 
12 1 1 and plotted vs the inverse system size L _1 in a log-log scale. Different symbols are used for different values of the power-law 
exponent a, as indicated on the plot. The case with only nearest neighbor interactions is also shown (N.N green crosses) for 
comparison. Right panel: Infinite size AF order parameter tuaf plotted vs a, obtained using finite size scaling of Stt(L) (shown 
on the left panel). The quantum phase transition between the Neel phase (mAF 7^ 0) and the QLRO phase (rriAF = 0) occurs 
at a c = 2.225 ± 0.025. The SW estimate (a" w ~ 2.46) is indicated by the arrow. 

to extract the thermodynamic limit behavior of S n , we perform a finite size analysis in order to get the AF order 
parameter, given by 



^SniL) -> m A F, (L— ► 00). 



(23) 
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Utilizing the fact that in the QLRO regime the spin-spin correlation functions decay as stated in Eq. Q , we therefore 
expect in this regime the following behavior for the staggered structure factor per site 

SAL) = \ j\-l) r (S ■ S r )dr ~ (ln ^ )3/ , (L-oo). (24) 

On the other hand, in the Neel phase, the finite size scaling of the order parameter can be evaluated using the small 
k SW spectrum (see appendix 0), leading to 

SAL)-m\ F ~L^ +0{L a - 3 ). (25) 

We used second order polynomial fits in L~ ~ to extrapolate the finite size data to their thermodynamic limit values, 
shown in the right panel of Fig. [21 The quantum phase transition between the AF Neel order and the QLRO phase 
is clearly visible for a critical value 2.2 < a c < 2.25. It is also interesting to compare this estimate with the one from 
SW approximation giving a™ ~ 2.46. 

Let us now concentrate on the a— and A— dependent Hamiltonian by keeping a fixed value for a while varying A. 
We first focus on the case with a = 2.1 which is expected to display a transition for a non-zero value of A. As pointed 
out by Reger and Young studying finite size AF clusters in d = 2 [2j|, the sublattice (infinite size) magnetization 
can be obtained either from the staggered structure factor [Eq. or from the correlation functions at the largest 
separation 

C(L) = (Si ■ S l+L/2 ) - ±m| F , (L->oo). (26) 

In the Neel phase, both estimators S n (L) and C(L) are expected to converge to m\ F with a similar power-law behavior 
but with different pre-factors. This feature is illustrated by the computation of C(L) and S V (L) for a = 2.1 and four 
different value of the long range term strength A = 3, 2 , 0.9, 0.6, as shown in Fig. [3J Using second order polynomial 



■ S,(L) ■ C(L) 




jja-3)/2 j^(a-3)/2 

FIG. 3: Staggered structure factor per site S W (L) [Eq. 1231 1 and mid-chain correlation function C(L) [Eq. 126H 1 computed at 
T — with QMC for a — 2.1 and 4 different values of A, as shown on the plot (a)-(d). Finite size scaling has been performed 
on finite systems up to L max = 128 for the largest values of A ((a) and (b)) and up to L max = 256 for the smallest ((c) and 

(d)). The dashed lines are polynomial fits of the form m\ ¥ + a\L~^~ + a2L a ~ 3 . 



fits in , we can obtain infinite size extrapolated values for ?tiaf(A) from S V {L) or C(L), as reported in table lTVAl 
Finite size effects are more pronounced for the staggered structure factor than for the mid-chain correlation function 
because S„(L) is the result of the integration of the staggered correlation function along the entire chain and therefore 
is sensitive to short distance terms. However, the estimates for the sublattice magnetization obtained from C(L) and 
S n (L) (see table lTVA"|l are both in good agreement, especially when the system is deeply in the Neel regime (large 
values of A). On the other hand, when the system is approaching the quantum critical point (QCP) where tuaf ^ 0, 
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A 


771AF from Stt(L) 


ttiaf from C(L) 


3 


0.353 


0.356 


2 


0.295 


0.301 


0.9 


0.106 


0.091 


0.6 





0.02 



TABLE I: Infinite size extrapolated values of the sublattice magnetization iuaf obtained for a = 2.1 and A = 3, 2, 0.9, 0.6 
from power-law fits of the staggered structure factor Stt(L) and the mid-chain correlation function C(L) (see Fig. [3Jl 



the finite size effects are significant enough to prevent us from obtaining a very precise estimate of the critical coupling 
A c where the AF LRO vanishes. 

Of course, in principle it is possible to perform very large scale numerical SSE simulations on the largest reachable 
system sizes, as we did for the A = 1 case with L max = 4096. However, since our goal here is to investigate the 
quantum critical phenomena in the A — a plane, we need a good sampling of this parameter space and we therefore 
restrict the simulations over systems of maximum size L max < 1024. We then use another strategy, based on scaling 
arguments, to perform a better data analysis close to criticality. This is described next. 



B. Scaling analysis 

As previously discussed, the finite size effects are bigger for the staggered structure factor than for the mid-chain 
correlation function. Therefore, we now focus on C(L) which is expected to saturate to a constant value in the Neel 
phase, whereas in the QLRO regime, the behavior C(L) — > y/\a{L/a)/L is expected, a being a non-universal constant. 
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FIG. 4: T = QMC results for the mid-chain correlation function C(L) computed for a = 2.1 and different values of A, as 
indicated on the plot, (a) C(L) is plotted vs the system size L for < A < 1.15. (b) Both x- and y-axis are rescaled using 
two parameters: the crossover length scale £a and ^a- The data collapse results in two universal curves: one for the Neel 
ordered phase (top one) and one for the QLRO regime (lower one). Note that for clarity, the QLRO universal curve has been 
shifted downwards. The red dashed-line materializes the critical separatrix between the two regimes, decaying with an exponent 
~ 0.63. Inset: Crossover length scale £a extracted from the data collapse in the QLRO (full circles) and Neel regime (open 
squares). The dashed lines are power-law fits of the form |A — A c |~" with A c ~ 0.45 (indicated by vertical dotted line) and 
v ~ 15. 



In order to illustrate the scaling analysis, let us continue to study the case with a = 2.1, as in the previous 
subsection. We have computed C(L) for several values of the long range coupling strength A in the range [0, 1.15], 
for sizes up to L — 256. The results, shown in Fig. 01(a), clearly show the existence of a finite critical value A c which 
separates the QLRO and the Neel regimes. In order to locate precisely this QCP, let us assume that a typical length 
scale £a governs a crossover from the QCP to the Neel phase if A > A c and from the QCP to the QLRO regime if 
A < A c . Precisely at the critical point, the spin-spin correlation function decay like a power-law 



C QCP (L) ~ 



(27) 
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thus defining the critical exponents rj and z, the critical dynamical exponent. Without making any assumption about 
the values of the aforementioned critical exponents, let us now define scaling functions f±(x), with x — L/^\, for 
A > A c and A < A c , respectively, by 



/±(s) 



C(L) 



(28) 



Hence, the scaling functions obey 



(29) 



X X 



l — Z — T] 



f_(x) cx .T- 2+2+r 'VhT^ if x > 1 and A < A c (QLRO) 
f+(x) cx x- 1+z+v if x > 1 and A > A c (NEEL) 
/±(0) = 1 if A - A c (QCP) 

It is convenient to also rescale the y-axis with the unknown function in order to get C(L) x vJ/ A — /(x) 
We then expect ^ A to be proportional to £ A +?,_ ■ Using such scaling forms, we have obtained the collapse of the 
data shown in Fig. 0Ja) into two universal curves shown in Fig. Efb) . The parameters £\ and *ff\ have been chosen 
to give the best data collapses. Using such a scaling analysis, we find a critical coupling A c = 0.45 ± 0.05 that we can 
compare to the overestimated value A c = 0.6 previously found using the more simple finite size scaling Eq. 1|25[) . The 
critical correlation (given by the separatrix between the two regimes in Fig.^b)) is characterized here by a power-law 
decay with an exponent (z + rj — 1)qcp — 0.63. Note also that the crossover le ngth scale £\, plotted in the inset of 
Fig. 0Jb), diverges on both sides of the transition with a large exponent v ~ 15 24] . These and other issues related 
to the critical exponents will be discussed in detail in Sec. IV CI 



V. QUANTUM MONTE CARLO II: PHASE DIAGRAM AND CRITICAL BEHAVIOR 



The scaling analysis described above has been repeated for several values of a in order to explore and construct the 
phase diagram of the model (@J in the A — a plane. 



A. a = 2: Marginal case 

Let us first focus on the marginal case with a = 2 for which a similar data collapse analysis is performed and shown 
in Fig. |S] for the mid-chain correlation function. 
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FIG. 5: T — QMC results at a = 2 for the mid-chain correlation functions C(L). As in Fig.0] both x- and y-axis have been 
rescaled in order to get the best data collapse. For the different values of A indicated on the plot, the data collapse on a unique 
crossover curve towards the Neel ordered phase. Inset: Crossover length scale £\ plotted in a linear-log scale vs A -0,4 . The 
dashed line is a fit of the form Eq. (gOJ with a = 0.4. 

Using QMC simulations results for chains up to L — 512 sites, with A € [0, 3], we have been able to get an universal 
curve (see Fig. [SJ) which shows a crossover towards a Neel order phase (i.e. C(L) — > constant if L ^> Note that 
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for A < 0.1, the typical length scale necessary to get a good collapse becomes very large so that there is no overlap 
between our different curves and the data collapse analysis is impossible to achieve. Nevertheless, the crossover length 
scale £a, plotted in the inset of Fig.0 displays an exponential divergence when A — > 0. Guided by the RG calculations 
presented in section l\*Tl we can fit the A-dependence of the crossover length scale by 

6,~exp(C7A CT ), (30) 

with o = 0.4 and C being a free parameter. It is however important to note that since £a suffers from large error 
bars, and so does the fitting parameter, we have forced a to its value found in Eq. (|66|l . 

Unambiguously, £a is found to diverge when A — * which means that at the marginal point a = 2, any A > will 
drive the system towards the Neel phase. In other words, the long range interaction perturbation of strength A is 
marginally relevant at a — 2. This result agrees with the RG calculations presented in section IVT1 



B. Phase Diagram 

As previously stated, when a < 2 the long range interaction is a relevant perturbation and any A > will drive 
the QLRO phase towards a AF ordered Neel phase with ttiaf 0. On the other hand, when a > 2 a simple power- 
counting tells us that the perturbation is irrelevant which should imply that the QLRO is stable against a small 
perturbation A > 0. It turns out that such a simple argument is not sufficient to provide a correct description of 




FIG. 6: T — phase diagram of the long range S — 1/2 model [Eq. computed by large scale QMC simulations and plotted 
in the A — a plane. A line of critical points (circles) separates a long range ordered phase (Neel) and a quasi long range ordered 
phase (QLRO). The error bars, due to some uncertainties in the finite size scaling analysis of the numerical data, are explicitly 
show on the plot. The dashed line is a guide for the eyes. 



the quantum critical behavior of the system (see the next section for a more advanced field theoretical description). 
Based on large scale numerical simulations, we provide hereafter a picture which is consistent with the existence of a 
non-trivial line of fixed points in the A — a plane. 

Using QMC simulations on systems of up to L = 1000, we performed the scaling analysis for the mid-chain 
correlation function C{L) as well as for the staggered susceptibility (see below) and computed the phase diagram for 
2 < a < 2.7. For each value of a, the QCP A c is found by the separatrix between the two crossover functions (see 
Fig. EJb)) with some error bars due to the discrete sampling in the A space as well as the strong divergence of the 
crossover length scale close to the critical point which makes the data collapse delicate. We present in Fig. [fj] the 
QMC phase diagram in the A — a plane. As discussed, A is marginally relevant at a — 2, driving the system towards 
a Neel phase with LRO. At small A, the critical line increases sharply from a = 2 and displays a negative curvature. 
By contrast, spin-wave theory (see Fig.^) and large- N approximation predict that a c (X — ► 0) = 1. In the range of A 
considered here (A < 8), the critical line stays well below the value a = 3 and we expect this feature to remain true 
for all A. This behavior is consistent with the proof of absence of LRO at T — for a > 3 [TEl ITt^. 
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C. Critical exponents 

The transition line between Neel LRO and QLRO is a non trivial line which displays continuously varying critical 
exponents, as we show now. 

1. Divergence of the crossover length scale 

The standard theory of quantum phase transitions involves a set of critical exponents which govern the universal 
behavior of various quantities close to or at the QCP. One of them is v which tells us how does the correlation length 
diverge in the real space direction close to the critical point. Usually this correlation length is defined in the disordered 
phase by the exponential decay ~ exp(— r/£) of the correlation function associated with the order parameter. In 
our case, the non ordered regime A < A c is already critical and thus the correlation length is intrinsically infinite. 
Nevertheless, the typical length scale £ A which governs the crossover phenomenon, diverges at the QCP (on both 
sides) with an exponent which we call v by analogy: 

koclA-Ad-". (31) 

As already discussed, an accurate numerical evaluation of the exponent v is difficult, because of some intrinsic 
uncertainties in the data collapse procedure. Nevertheless, at our level of precision we observe this crossover length 
scale exponent increasing when a — > 2 + . In particular, at the marginal point a = 2, we find an exponential divergence 
of £a near A = [Eq. (|30|) ]. formally corresponding to v — oo. This divergence of v when a approaches 2 is actually 
in good agreement with the results of field theory and RG calculations presented in Sec. IVII 

2. Staggered magnetization exponent and hyperscaling relation 

The other scaling parameter ^ \, used for the collapse of the correlation functions data, also contains some infor- 
mation. First, if the scaling hypothesis used above with the help of the crossover functions is correct, we expect 
* a ~ £a + " _1 which 

gives another estimate for the critical exponent of the decay of the correlation function [Eq. (|27l) ]. 
This is illustrated in Fig. [7| where is plotted vs £\ for a = 2.1. Data, presented for both sides of the transition 
in Fig. \7\ clearly display power-law dependences with an exponent r\ + z — 1 in very good agreement with the value 
of 0.63 previously found along the separatrix in Fig.|2fb). Note also that the agreement is even better when getting 
closer to the QCP. 

3 
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10° io' io 2 io 3 io 4 io 5 

FIG. 7: Log-log plot of the scaling parameters and £a obtained from the collapses shown in Fig. |ljb) for a — 2.1 in both 
phases: Neel (O) and QLRO (□). Within error bars (explicitly shown on the plot), data are fitted by power-laws (dotted lines) 
of the form $ 617 for the Neel regime and ^J' 636 for the QLRO. 

In the ordered phase, according to Eq. J2EJ), we expect for the AF order parameter 

m-AF oc £ A 2 cx (A - A c ) 2 . (32) 
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This implies the usual hyperscaling relation involving the critical exponent (3 governing the onset of the order parameter 

m A F oc (A - \ c f, (33) 

which must therefore satisfy: 

2j3 = v{z + rj - 1), (34) 

i.e. the usual hyperscaling relation in d = 1. 

3. Analytical estimate of the exponent r\ 

Following the same philosophy as in the mean field argument given in section [HI we can calculate the expectation 
value of the long range perturbation 

J(r,r')(S r ■ S r ,) ~ lr _ r , ] l +z+n _ 1 , (35) 

at the QCP, with some unknown critical exponents z and rj. The finite size correction to the free energy density now 
scales like £ 2 ~ Q ~ Z ~' ) . The singular part of the free energy at some non trivial QCP is expected to scale like L~ 1_z 
for a finite size system. Thus the two corrections will scale in a similar way if 

X] = 3 — a. (36) 

The same condition is also obtained by demanding that the long range interaction be invariant under the RG trans- 
formation involving a scale factor s: 

n(r, x) S - {z ~ 1+rl)/2 n{T/s z , x/s). (37) 

(The rescaling factor of s _ ( z_1 +'')/ 2 for n implies the equal time correlation exponent of z — 1 + rj.) Rescaling x and 
r inside the integral Eq. (|15fl (which represents the contribution of the long range term in the action), the condition 
for invariance under this RG transformation leads to 2 + z — a — (z — 1 + r/) =0, which gives the expression (|36|l for 
rj. As already shown in section ITTT1 Large- N calculations also give the same value for r\ [2f|. Let us mention that the 
RG analysis, presented below in section IVll also agree with such an estimate, up to order (a — 2) 2 . 



4- Numerical determination of the exponent r\: scaling of the staggered susceptibility 

The T = staggered susceptibility, defined on a finite ring of size L by 

1 f°° 

XW- Z EH) NI J (SM-S^r^dr, (38) 

ij 

obeys the standard finite size scaling at the QCP: 

x(tt) cx (39) 

Also we know, for instance from SW calculation (see appendix [XJ, that in a Neel ordered state the staggered suscep- 
tibility will scale quadratically with the size L. On the other hand, in the QLRO characterized by r\ = 1, we rather 
expect a linear scaling of x(7r) with L. Consequently, there are three distinct regimes for the staggered susceptibility: 

{constant if A > A c (NEEL) 
L~ x if A < A c (QLRO) (40) 
if A = A c (QCP). 

We use the same scaling procedure as for the correlation functions, to obtain data collapses onto two different curves, 
as illustrated in Fig. |H1 for a — 2.2. Indeed x(7r), computed with QMC on chains of up to L = 1000 sites, displays 
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FIG. 8: QMC results for the T = staggered susceptibility x( 7r ) [Eq. jSBj ] computed for a = 2.2 on systems up to L = 1024 
spins. The upper inset shows x( 7r ) vs L f° r various A £ [0.1, 1.9]. The main plot shows the results of a data collapse onto two 
universal curves, after a rescaling of both x- and t/-axis using two parameters £a and 0\. Asymptotically, the LRO curve (top 
one with data for 0.95 < A < 1.9) saturates towards a constant whereas the QLRO one (lower one with data for 0.1 < A < 0.85) 
displays a L _1 behavior, characteristic of r\ = 1. Between them the separatrix shows the critical behavior around the transition 
at A c ~ 0.9, decaying like £ -0 - 8±001 . The lower inset shows a log-log plot of the scaling parametrs Q\ and £\ used to achieve 
the collapses in both phases, Neel (O) and QLRO (□). Within error bars, data can be fitted for the entire range by power-laws 
(dotted lines) of the form £° x 785 for the Neel regime and 815 for the QLRO. 




FIG. 9: Numerical estimate for the critical exponent r\ along the critical line obtained using the critical behavior of 
computed with QMC (red symbols). The red dashed line is a guide to the eyes. The blue full line is the analytical estimate 
r\ — 3 — a. 



clearly a crossover phenomenon on both sides of the transition, also characterized by a crossover length scale which 
is directly proportional to the one previously extracted in the analysis of the correlation functions. Note also that 
such an analysis provides a second physical observable way to locate the QCP: in fact, the analysis of C(L) and 
x(tt) both agree (within the error bars) on the value of A c . Moreover, we expect the scaling hypothesis to be valid if 
®\ ~ This is actually the case, as illustrated in the lower inset of Fig. [H] where we find i] = 0.8 ± 0.015. We can 
also obtain the quantum critical exponent r\ from the separatrix between the Neel and QLRO regimes (see Fig. El, 
which is expected to decay as L~~ v . For a = 2.2, we find 77 = 0.8 ± 0.01. We have repeated this computation of 
x(tt) for several other values of a £ [2.1,2.7] to calculate the corresponding r](a). The results are plotted in Fig. El 
and compared to the previously discussed estimate 77 = 3 — a. It is very remarkable to see how this rough estimate 
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reproduces quite well the actual value. Only for a > 2.3 a deviation starts to appear. 



5. Dynamical exponent z < 1 

The dynamical critical exponent z, involved for instance in the critical decay of the correlation function Eq. J37J, 
can be evaluated from the spin-spin correlation function at the QCP. From the fit of the separatrix between the two 
data collapses (see for instance Fig. Q]where for a = 2.1 we estimated r\ + z — 1 = 0.63 ± 0.03) we obtain an estimate 
for rj + z — 1. Then, using the estimates of 77, determined separately with the staggered susceptibility, we obtain a 
numerical evaluation of z. Results are shown in Fig. ^| for 2 < a < 2.7. For a = 2, the QCP at A = displays the 
critical behavior of the short range model, with 77 = z = 1. Surprisingly, when moving from a — 2 along the transition 
line, z becomes very rapidly < 1 and, within the error bars, seems to saturate around a value ~ 0.75. It is actually 
natural to expect z 7^ 1 since the long range interaction breaks Lorentz invariance. However, unlike for the estimate 
of 77, the dynamical exponent obtained within the large- iV expansion z = (a — l)/2 does not agree with the QMC 
results. As we discuss in the next section, using a "RG improved" perturbation theory, z is found to be < 1 but does 
not display such a big reduction. 

^ — ■ — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — ■ — 1 — ■ — 1 — ■ — 
I 

I 
\ 



"' J 2 ' XI ' 2^2 23 ' 2.4 25 2.6 2.1 ' 

a 

FIG. 10: Numerical estimates for the dynamical critical exponent z along the critical line obtained using the critical behavior 
of the correlation function [Eq. 1271 1 computed with QMC, and the numerical estimate of n (see Fig.|5Jl. The numerical results 
(open circles) suffer from large error bars, as shown on the plot. The green dashed line is a guide to the eyes. 
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VI. FIELD THEORY /RENORM ALIZ ATION GROUP RESULTS 

The low energy, continuum limit, imaginary time action takes the form: 

S[n] = S [n] - g [ drdx^r, x) - \a a ~ 2 f ^^H(r, x) ■ n(r, y). (41) 
J J \x-y\ a 

n(r, x) is the antiferromagnetic order parameter field defined by the continuum limit expression of Eq. I|13f) . Here 
So is the action for a free massless relativistic boson, in terms of which n may be represented in a non-linear way. 
Equivalently, we may regard it as the action of the k — 1 Wess-Zumino-Witten non-linear a model. The field <j> is 
defined as: 

= 2ttJ l • J R , (42) 

and is normalized so: 

< 0(*)0(O) >= t^— JTTJ ■ (43) 



The corresponding coupling constant, g has a bare value of order 1 for the short range AF chain and is marginally 
irrelevant. It is responsible for various logarithmic corrections such as the one in the correlation function of Eq. |J3}. 
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Note that the dimensionless coupling constant for the long range interaction A is only proportional to the one used 
before for the lattice microscopic model in Eq. Q, and a is a short distance cut-off with dimensions of length. As 
already noted in section^ since n has a scaling dimension of 1/2 from Eq. 114JI . A is irrelevant for a > 2, relevant for 
a < 2, and marginal for a = 2. Also note that A > corresponds to non- frustrating interactions which favor the Neel 
state where (n z ) ^ 0. Our strategy is to study this model when < a - 2 « 1 and < A <C 1 using perturbation 
theory in g and A. Since g renormalizes to at large length scales, when A = (i.e. in the short range model) this 
can give useful results, for small bare A, even when the bare value of g is 0(1). We will find that an interplay between 
the local marginal coupling constant g and the irrelevant non-local coupling constant A governs the critical behavior 
in this regime. 

We now consider the low energy effective field theory for the long range model in Eq. (|41|l . in the limit of small g 
and A, using RG methods. When A = 0, the RG equations reduce to the standard ones for the short range model. 
These take the form: 

=-. 9 2 -(l/2).g 3 + .... (44) 



din a 



Here we define our RG transformation by increasing the short distance cut-off a. The bare value of g is positive for 
any non- frustrated short range model and is typically 0(1). The basin of attraction of the g = fixed point is known 
to extend to such large bare values of g so that g = is the universal stable fixed point for short range models. The 
flow of g towards zero at long length scales is controlled by the quadratic term in the /?- function of Eq. I|44|) . giving: 

9(a) -> w v (45) 
ln(a/a ) 

where ao is the original cut-off and a is a larger value obtained from integrating out modes with wave-lengths between 
ao and a. This logarithmically slow flow of g(a) to zero is responsible for logarithmic corrections to the correlation 
function and other properties of the short-range models. A linear term in the /3-function for A follows immediately 
from the factor of a a ~ 2 in Eq. 1)41(1 which in turn is a consequence of the fact that n has scaling dimension 1/2: 

' /A (2 - a)X + . . . (46) 



din a 



So, ignoring the effects of g, A grows larger at long length scales for a < 2 but smaller for a > 2. Long range 
interactions are irrelevant for a > 2. However, it is necessary to consider higher order terms in the /3-functions for 
both g and A to understand the phase diagram, even at a « 2. 

To calculate additional terms in the /3-functions, we define our ultra-violet cut-off by forbidding any 2 points in 
space-imaginary time from getting closer than a, in a perturbative calculation of the partition function (or long 
distance Green's functions). In particular, this means that the long range term in the action is cut-off as: 

S X = -Aa- 2 / ^|n(r, x) ■ n(r, y). (47) 

J\x— y\>a !•£ V\ 

When the cut-off is increased from ao to a — ao + 5a, there is an additional change in S of first order in 5a, which 
comes from the change in the integration region: 



5S = — Xa 



a-1 





' p—ao 




j dxdr 


L + 


n 






Ja J 



-. — ^—nir, x) ■ n(r, x + u). (48) 
\u\ a v ; v ; K ' 



Since both factors of n are very close together, we may use the operator product expansion. This follows from the 
3-point Green's function: 

< n a ( Zl )ct>(z 2 )n b (z a ) >= ~ iggl (49) 
This implies the operator product expansion (OPE): 

n a (z)n & (O)->^|z|0(O) + ... (50) 

Using this in Eq. (|48l) . gives: 

5S = -Xa*- 2 j d T dx(t)(T,x)AT: f — ^ sa -^4ttA j drdx^r, x). (51) 
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This corresponds to a renormalization of g: 



Sg = 4-7T — A 
a 



and hence to another term in the /3- function for 5: 



din a 



= 4ttA - .g 2 



(52) 



(53) 



There is one more term in the RG equations that is important at small a — 2, corresponding to a correction to A of 
order Xg. This can be calculated from the OPE: 



(j)(z)n a {0) 



1 



8tt|;; 



rn a (0) 



giving: 

dX/dhxa = -{a - 2)A + {l/2)Xg + ... 
The RG equations, Eq. H53[l and (|55|l have an unstable fixed point for a > 2, at: 

g c ffc I'm - 2) 
A, 



(54) 



(55) 



2(a 

±(a-2) 2 . 



(56) 



For a < 2, a positive A always runs away to large values as we lower the cut off (i.e. increase a), corresponding to 
LRO. On the other hand, for a > 2, a small enough positive bare A flows to zero while a larger bare value flows to 
large values (see Fig.^Jl. These statements remain true even when the bare value of g is 0(1) as we expect it to be 




FIG. 11: Renormalization group flow of Eqs. 153H and 155H in the case a = 2.3. The dotted line represents, schematically, the 
values of the bare couplings in the field theory as the parameter lambda in the lattice model is varied. The unstable fixed point 
at A c = 0.0286, g c = 0.6 separate the flow to the stable fixed point at A = g — which represents the quasi-long-range ordered 
phase and the flow to infinite A, g, which represents the long range ordered phase. The black lines with double arrows denote 
the separatrixes between these two phases. (The corrections to the flow equations are presumably significant for this large a 
value of a — 2, but we graph this case for ease of visualization. 



in general for a short range spin chain. For a small bare A, g initially renormalizes towards small values as it would in 
the short range chain until eventually Eqs. (|53|l . (|55(l becomes valid. The stable A = g = fixed points corresponds 
to the standard QLRO phase of the short range spin chain. The non-trivial unstable fixed point separates the ordered 
and quasi long range ordered phases. Of course there are higher order terms in both RG equations, but they do not 
invalidate our conclusions on the location of the fixed point, for small enough a — 2. Both terms on the right hand 
side of Eq. (|53|) are 0[(a — 2) 2 ] at the fixed point; any possible higher order terms such as g 3 or A 4 are at least of 
0[(a — 2) 3 ]. Similarly, both terms on the right hand side of Eq. (|55l) are 0[(a — 2) 3 ] at the fixed point; higher order 
terms are at least of 0[(a — 2) 4 ]. To reach this conclusion it is important to realize that there cannot be any terms in 
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dX/dhia which contain no factors of A; a purely short range interaction cannot generate a long range one although 
the reverse is not true. 

Thus we appear to have a rare example of a non-trivial fixed point which can be accessed perturbatively. (But 
see the discussion below of potential problems with this app roach.) We note that a similar expansion for long range 
classical spin models was introduced in 0. See also |2?l l28j. Our quantum spin model, in the continuum limit, non- 
linear a model approximation, appears rather similar, in the imaginary time path integral formulation. An important 
difference, however, is that our model has an action which is long range in the space direction but short-range in the 
time direction. Thus it corresponds to a classical model in 2 space dimensions with short range interactions in one 
direction and long range interactions in the other. It is this asymmetry which leads to a dynamical critical exponent 
z < 1. Another important difference from a 2-dimensional Heisenberg model is the topological term in the short- 
range part of the action which is responsible for the quasi-long-range order. We remark that an integer-spin quantum 
Heisenberg chain with long range interactions could be expected to have identical critical behavior to a classical 
Heisenberg model in two dimensions with interactions which are long range in one dimension. We also remark that 
an xxz quantum spin chain with long range interactions could be expected to have the same critical behavior as a 
two-dimensional classical xy model with interactions which are long range in one dimension. 

The phase boundary (or separatrix) can be found by determining the line in the g — A plane which renormalizes to 
the critical point. Combining Eqs. (|53|l and l|55(l gives: 

K dX . . . f 90 dg(g-g c ) 

„ ~ = {1/2) L 9 ^9Kx/x c y (57) 

where A = X(g) is a function of g along the RG flow in the integral on the right hand side of Eq. I|57[l . go and Ao are 
the values at arbitrary points on the separatrix. Since A increases monotonically to the value A c with increasing a, 
we may obtain an upper and lower bound on the right hand side by replacing X(g) by A c and respectively inside 
the integral on the right hand side of Eq. I|57|) : 

(1/2) ln(.g /e.g c ) + g c /2g < ln(A c /A ) < (1/2) ln[(g + g c )/2g c ] . (58) 

Now using the fact that g c << go, this becomes: 



Vgo/eg c < Xc/Xo < V W 2 ffc, (59) 

that is: 

.637(a - 2) 5 / 2 /^ < < -743(a - 2) 5 / 2 /^. (60) 

Thus on the separatrix Ao is O[(o, — 2) 5 / 2 ]. [A numerical solution of the RG equations indicates that Ao is very close 
to the lower bound in Eq. JBDJ.] Assuming a bare g of 0(1), it is then possible to predict the shape of the phase 
boundary in the A — a plane close to a = 2. There is an unknown multiplicative factor relating the lattice coupling 
A to the continuum coupling A. However, we can predict that 

a c {X) -» 2 + CA 2 / 5 , (61) 

for some unknown constant factor, C, as a — > 2. As mentioned above, we expect that 2 < a c (A) < 3 for all A and all 
S. 

While our QMC results also predict that a c — ► 2 when A — * 0, but otherwise Eq. H61fl does not agree well with the 
QMC result as shown in fig. ^] It is interesting to note that lowest order SW theory and the large- N approximation 
make the mean-field prediction that a c — > 1 as A — > 0, in clear disagreement with our RG and QMC results. 

Linearizing the /3-functions at the critical point gives: 




(a-2) 2 /(27r) \ / A - A c 
4tt -4(q - 2) )\g-g c 



(62) 



This matrix has one positive (unstable) right eigenvalue, (\/6 — 2) (a — 2), implying a crossover length scale: 

^oc\X-X c \- v , (63) 

with a critical exponent: 

1 



(VQ-2)(a-2Y 



(64) 
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FIG. 12: QMC phase diagram at small A compared to the RG prediction of Eq. 1)61^ with C = 0.14 and C — 0.11. 



which diverges as a — > 2. The corresponding unstable direction is 

(a -2) 



A — A n = 



2tt(V6- 2; 



(9- 9c)- 



(65) 



For a = 2 exactly, LRO occurs for any A > but behavior characteristic of the quasi long range ordered fixed point 
occurs out to a cross over length scale: 



£ cx exp[C/A 



2/51 



(66) 



for a constant factor C. 

We may also determine the critical exponent, 77 + z — 1, controlling the equal-time correlation function at the 
non-trivial critical exponent: 



(67) 



when a — 2 << 1. Since A and g are both small at the critical point, for a — 2 << 1, we may simply do "RG improved" 
perturbation theory. That is we calculate the correlation function to first order perturbation theory in g, replace g 
by g c , and interpret the result as the expression Eq. I|67(l with 



T 1 + z = 2 + 0(g c ). 



(68) 



[To lowest order in (a — 2) we only need consider g, not A, since A c cx g 2 c \ The Green's function, up to first order in 
perturbation theory in g, the bare coupling, is: 



\z\ 2y/3 



= ^-r + -^8 ab I d 2 z' 



z' r z - z' 



(69) 



The integral in Eq. (|69(l must be restricted to the region \z'\ > a, \z — z'\ > a. For \z\ >> a, the integral is dominated 
by the two regions, \z'\ « \z\ and \z' — z\ « \z\, giving 



< n a {z)n b {0) > m — 
\ z \ 

gab 

Now replacing g by g c , its fixed point value we obtain: 

g ab 



2^3 J \z'\ 



l + -^=47rln(|2|/a) 



2V3 



< n a (z)n b (0) >= — {1 + (a - 2) \n(\z\/a) + 0[(a - 2) 2 ]} 



(70) 



(71) 
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Now, using the fact that the correlation function should have a pure power-law form at the fixed point, we may 
interpret this result as the leading term in the expansion of: 

< n«(z)n\O)>= J -0^. (72) 

To this order in (a — 2) we obtain the same exponent for z = ix or z = r, implying that the corrections to the 
dynamical exponent, z = 1, are higher order in (a — 2). Thus: 

V = 1 - {a - 2) + 0[(a - 2) 2 ]. (73) 

This is the same value of r\ found in the large N approximation and by a simple scaling argument in sub-section 
(V-C-3). It agrees well with QMC results for a < 2.3 as shown in fig. EI It is natural to expect that z ^ 1 since the 
long range interaction breaks Lorentz invariance. Since u> oc |fc|^ a_1 ^' 2 in the ordered phase it is natural to expect 
that z < 1 at the critical point, as is found in the large N approximation. However, clearly 1 — z must be at least of 
0[(a — 2) 2 ] since the long range coupling constant, A is of that order. In fact, we suspect that 1 — z is even higher order 
than quadratic in (a — 2). This conclusion does not fit well with the QMC results for z, presented in fig. El There it 
was found (although with large error bars) that z appears to have a nearly constant value, z « .75, for a > 2.1. As 
a is further decreased z appears to rise very rapidly towards 1. 
So far, we have ignored another possible interaction: 

S- (s'tt/3) / dTdx{Jl + J%). (74) 



This interaction is, in fact, present for the short-range spin-chain with a large coefficient. Since the Hamiltonian 
for the k = 1 WZW model can be written quadratic in currents, this "interaction" term can be regarded as simply 
shifting the velocity, which we have so far set equal to 1, to: 

v^l-g'/2. (75) 

The RG equations, for the short range model, including g' take the form, to cubic order: 

dg 



din a 
dg' 
din a 



-3 2 -(l/2) 3 ( 5 2 + 5 ' 2 ) (76) 
(3/4) ff 3 . (77) 



Starting with g, g' > and 0(1), these equations predict that g — ► and g' flows to a value of 0(1). The large value 
of g' at the fixed point can simply be interpreted as a large renormalization of the velocity, provided that g < 2. 
In fact, this is what happens, for example in the Hubbard model at half-filling. The spin velocity is reduced by the 
Hubbard interactions. An alternative approach is to adjust v to the correct value and drop g' completely from the RG 
equations. In fact, Eq. (|77fl depends strongly on the renormalization and cut-off scheme. With a Lorentz invariant 
cut-off and renormalization procedure, the non Lorentz-invariant term, proportional to g' will not be generated under 
the RG if it is initially absent. Breaking of Lorentz invariance in this problem at low energies, just means shifting 
the velocity. If we work directly with the exact velocity, then it is apparently permissible to set g' — and use 
a Lorentz invariant renormalization procedure so that g' remains zero under renormalization. In fact, using this 
procedure in Eq. (|76|) leads to various predictions of logarithmic corrections which are in good agreement with Bcthe 
ansatz and numerical results for the S=l/2 chain. If we set g' equal to some arbitrary non-zero value in Eq. I|7t)|) the 
coefficients and powers of log corrections would change, resulting in worse agreements with numerical results. Thus, 
this procedure of ignoring g' seems to be a valid and useful one. 

We now consider the interplay of the long range coupling constant, A, with the non-Lorentz invariant local coupling, 
g' . The needed OPE's can be obtained from the general conformal field theory result for the 3-point Green's function 
of the energy momentum operator with a primary field of left-dimension 1/4: 



<T(zK(z 1 )n b (z 2 )>^^^ 

»=l 



1/4 1 d 



(z — Zi) 2 z — Zidzi 



Sjab 

(78) 



21 - Z 2 \ 



Here T = (2?r/3) J£, is the left -moving part of the Hamiltonian. This gives: 



< T{z)n a { Zl )n\z 2 ) >= ~ ^ gf rj. (79) 

87T \zi - z 2 \{z- z 1 ) I {z - z 2 ) 2 



19 



Also using the 2-point function of T: 



we can deduce the OPE: 



<T(z l) r(, 2 ))> = 2(27r)2(zi _ z2)4 ; (80) 



n a (z)n b (0) - ^pT(O). (81) 

Now consider the case where the separation is in the space direction, z — ix 

n a (x)n b {0) -> 5 ab n\x\[(2/3)(j)(0) - T(0) - T(0)] + .... (82) 

Here we have included the identical OPE coefficient for T = (27r/3) J^. and also the coefficient deduced earlier, in Eq. 
(|50|l . Using the same cut-off and RG transformation procedure as above, this implies a term in the /3-function: 

^f- = -6nX. (83) 
am a 

This drives g' towards negative values, corresponding to increasing the velocity. If we continue to use our previous 
RG transformation, we find that the dA/dlna does not pick up an term oc Xg'. The difference from the non-zero Xg 
term in Eq. (|86[) arises from the fact that the OPE is now: 

T(zK(0) - rtAO). (84) 

8lTZ Z 

The RG transformation gives: 



<5Aoc J — , (85) 

where, as before, the integral is over a circular shell with radius between a and a + 5a. This integral vanishes by 
rotational invariance. Thus the complete set of RG equations to low order is: 



dX 
din a 

dg 
din a 

dg' 
din a 

These equations have an unstable fixed point at: 



-(a-2)X + Xg/2 

4ttA - g 2 - (g + g')g 2 /2 

-6ttA. (86) 



A c = 

g c = 2(a-2) 

g' c = -2-2(a-2)w-2. (87) 

In this approximation, we obtain the same prediction for r] w 1 — (a — 2), as before and still get z ~ 1. There is a 
shift in the velocity of 0(1). 

The large value of g' at the fixed point makes the predictions of this RG analysis intrinsically suspect. As in the 
short-range case, we might agree to set v equal to its renormalizcd value and then drop g' from the RG equations. 
This leads to the same predictions about the value of A on the separatrix, and v as obtained above. 

However, there are some worrisome features of this RG analysis which arise from the long range interaction. We 
expect a dynamical exponent z < 1 at the critical point. It then does not make sense to use a rotationally invariant 
(i.e. Lorentz invariant) RG transformation. We would then get back a g 3 term in dg' /din a and, very importantly, a 
Xg' term in dA/dlna. We would then generally find that g is not small, 0(a — 2), at the fixed point. A would also 
not be small at the fixed point. In this case we would lose all perturbative control over the critical behavior even for 
a only slightly larger than 2. In this case, the unstable critical point would not be close to the QLRO fixed point, for 
a close to 2. One possibility is that the effects associated with z < 1 can be ignored to lowest non-trivial order a — 2. 



20 



Then our use of a Lorentz invariant RG transformation may be justified, fn this case, the fixed point really is close 
to the QLRO critical point for a sightly larger than 2 and our predictions for rj and v are correct in this limit. 

The QMC results seem to give at least partial confirmation of the validity of an RG approach based on the a — 2 
expansion. Most importantly, a c — > 2 as A — > and the critical exponents ij and z appear to approach their values 
in the quasi long range ordered phase (77 = z = 1 ) in this limit, with v diverging. Furthermore, excellent agreement 
with the prediction for 77 was obtained, over a rather large range of a (up to 2.2), as shown in fig. El We were not 
able to obtain accurate estimates of v from QMC to test the RG prediction. On the other hand, z showed rather 
surprising behavior, in fig. 1101 dropping rapidly from 1 to about .75 as a is increased from 2 to 2.1. This suggests 
that the asympotic, small a — 2 behavior may only occur for very small values of a — 2 << .1. Numerical difficulties 
preclude obtaining QMC data in this region. Furthermore, the phase boundary as determined by QMC, a c (X) could 
not be fit well to the RG prediction of Eq. Ijtilfl . except possibly at very small a c (where we have no data) as seen 
in fig. 1121 This could be interpreted as meaning that our RG approach based on an a — 2 expansion is correct in 
principle but is only valid in practice for extremely small values of a — 2. (The good agreement for rj then appears 
fortuitous.) Alternatively, the discrepancies may indicate a problem with our RG approach, perhaps resulting from 
our cavalier treatment of the non-Lorentz invariant interaction in Eq. (|74|l . 



VII. CONCLUSIONS 



We have studied long range non-frustrating 5 = 1/2 antiferromagnetic chains using spin- wave theory, large- N 
approximation, quantum Monte Carlo and analytic renormalization group methods based on an expansion in a — 2. 
All methods predict to a line of critical points in the A-q plane with continuously varying critical exponents. This 
critical line separates phases with true Neel long range order and quasi long range order. Quantum Monte Carlo and 
renormalization group methods indicate that this critical line terminates at A = 0, a = 2 and suggest that, along the 
critical line as a — > 2 + , rj w 3 — a, while v diverges and z — > 1 _ . 
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APPENDIX A: CALCULATION OF FINITE SIZE CORRECTIONS FROM SW: CONTRIBUTIONS 

FROM THE k = AND FINITE k MODES. 

1. General method 

Let H(h) = H — hO, where O is an operator and ft. is a field. If we denote by \h) the GS of H(h), it follows that the 
GS energy is Eas{h) = (h\H — hO\h). Since (h\h) = 1, it is straightforward to show (in direct analogy to Feynman's 
theorem) that: 

(d) = {h\d\h) = -™%@. (Ai) 

In general we need the expectation values of various operators O in the unperturbed GS, i.e. in the limit h — ► 0. It 
follows that all we have to do is to compute the change in the ground-state energy, due to the perturbation —hO, to 
first order in h. 

In the remainder of this Appendix, the Hamiltonian Ti is that of Eq. (0J. We are interested in finite-size chains 
with an even number of sites L, and periodic boundary conditions. 
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2. Staggered susceptibility 



Let O — J2i(~^) l Si ■ Then, according to Eq. ijAlf) . the staggered magnetization at T = is: 

m w = E(-i)^> = - 



/)— o 



and therefore the staggered susceptibility is: 

xO) = 



1 dM„ 



1 d 2 E GS 



L dh 2 



a. TTie fe = contribution 



We Fourier transform the spin operators, Sa n = 2/L^ fc exp {ik(2na))S^ 1 San+i = 2/L^ fc exp (ik(2n + l)a))S£, 
and collect only the fc = components. Let us denote Si = Sj: =0 = ^ n S2 n and S2 = S£ =0 = S-zn+i the total 
spins of the two magnetic sublattices, of L/2 spins each. Since we are in a Neel ordered state, Si and S2 are spins of 
total magnitude Sl — LS/2 = L/A for spins S = 1/2. Then, up to some constants that do not depend on h: 



where 



Hk=o(h) = Hk=o — Hi 



Ttk=o = jSi • S2, 



(A2) 



(A3) 



and 



L 



Hi =/i(Si 2 -S|) 



1 + \Y- 1 —— 

n>l y ' 



2 J, 



eff 



(A4) 



(A5) 



We need to do perturbation theory to second order in h, to find the staggered susceptibility. The ground-state of 
Hk=o is the state: 



\0) = \St = Q,M t = 0,S l ,S l ) = -=±= jr (-l) m \m,-m) 

V 7U— — Sl 



(A6) 



where St = Si + S2 . The perturbation links this only to other states with Mt = (see below) . Let us denote 

\S T ) = \S T ,0,S L ,S L ) (A7) 
where St = 0, 1, 2S^. With this notation, the second-order correction to the ground-state energy is: 

2 



However, 



2h 



Sl 



1 1 V2S L + 1 
Interestingly enough, one can show that: 



(A8) 



(A9) 



|l) = ai Yl (-l) m m\m,-m) 

m=—St 



(A10) 
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where the normalization constant is: 



ai = p L (s L + i)(2s L + iy (A11) 

It follows: 

(l|«i|0) = (A12) 
V2dl + i-a.1 

and (n|Wi|0) = 0, Vn > 2. By direct calculation, we find: 

E = -jS L {S L + 1); E 1= j- jS L {S L + I) (AI3) 

and therefore: 



(2) _ \(l\Hi\0)\ 2 M 4h> 

A ^ GS " "^b^T - -(25 L + l)a?j - ~^ ( ^ + 1} (AU) 



Since j = 2J eff /L, S L = L5/2, we find: 



As a result, the contribution of the k — modes to the staggered susceptibility is: 



1 d z E GS 



L dh 2 



2 (L 



h=0 ^e// 



3 Jeff V 2 



—5 + 1 ~ Zr (A16) 



b. The contribution of finite k modes 

Within the spin-wave approximation, the contribution of the k ^ modes to the ground-state energy of TL — hO is 
the zero-mode energy: 

E GS ~ ]T io k = £ V(l-f(k) + hy-g 2 (k) (A17) 

Here, the spin-wave dispersion is changed by the addition of the perturbation — h ^Y^i- The second derivative 
of Eqs with respect to h can now be calculated trivially, and in the limit h — > we find: 

Li — UJ, 

fc k 

where g(k) was defined before Eq. (|17fl . As A; — > 0, #(fc) — * 7 = const, u)k ~ &~ 2~ , and therefore 

Xfe#o(7r)~- 2 (A19) 

If a < 3, this is a smaller power than the L 2 contribution obtained from the k = mode. It follows that within the 
Neel ordered state, the staggered susceptibility scales like 1? (at least within the spin- wave approximation). 

3. Transverse correlation function 

We now choose 

= J2 S tnSr (A20) 
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so that its expectation values (6) is the transverse contribution to C(n). For simplicity, we assume n to be even 
(calculations can be done similarly for odd n). Since this is not a Hermitian operator, let: 



(A21) 



(A22) 



Both these operators are hermitian. According to Eq. (|A1|) and using invariance to translations, we then have 



L(S+S„) = (J2 S+ +n Sr) = (6 a) - i(0 B ) -> (S+So) =-j 



dE GS i dJ ^GS_ 



dh 



dh 



(A23) 



where E^J B ' are the ground-states energies in the presence of perturbations —hC>A/B- 



a. The contribution of finite k modes 

After a Fourier transform, we use the Holstein-Primakoff representation for all k ^= modes. Keeping only quadratic 
terms, we find: 



6a = 4S*^ cos(nfca) (b\b k + a\ 



i k a k 



(A24) 



B = sin(nfca) [b\b k - a\a k 



(A25) 



After adding this to the unperturbed Hamiltonian (in the SW approximation) and diagonalizing, we find the 
ground-state energies to be: 



E< GS = ^k.A + 4s /i cos(ro£;a) 

k=£0 fc^O 



(A26) 



EqJ = JS y ^ LUk.B — 4s/i y ^ sm(nka) 



■'E< 



i E> 



where 



.1 = \/ [7 - f(k) - y cos(nfca)] 2 - [g{k)] 2 



(A27) 



(A28) 



"k,B = \/h- f(k) + — sin(nfca)]2 - [g(fc)]2 



(A29) 



After taking the first derivatives and setting h = 0, we find the k 7^ modes' contribution to the transverse 
correlation to be: 



(q+q -,_2S^ (7 - f(k))e* k ™ 2S 



25' x - 

— E' 



(A30) 



In the limit fc — > 0, /(£;) — > 0, — > fc 2 and therefore: 



(S+S ~) = ^^(1 + ...) + 2S/L 



(A31) 



The second term is the second sum ( X^fe^o e " = $n,oL/2 — 1 = —1, since n > 0). For a > 2, the L 2 term is 
dominant. 
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b. The k — mode 



Keeping the full k = contributions, we find 



2h 



n k=0 (h) = jSi -s 2 - T {s+Si + s+Sz 



(A32) 



The notation has been introduced in the previous section. The ground-state |0) of Tlk=o is known (see section on 
staggered susceptibility), so the first order contribution to Ecs(h) can be evaluated directly: 



L(S^S ) — 



The calculation is trivial, and we find: 



dE, 



GS 



dh 



h=0 



(A33) 



2S 2 



2 

Is 



1 

L 



(A34) 



It follows that for this correlation function, the finite k modes give the dominant L dependence, which isL 2 . 

This calculation can be repeated for the parallel contribution to the correlation, (S*Sq). The L dependence remains 
the same, so we conclude that in the Neel state and within SW approximation, C(L) ~ . 



APPENDIX B: LARGE-TV CALCULATION 

Considering the order parameter <fi is a iV-component unit vector field 

\$(r,x)\ 3 = l, (Bl) 

the action can be written like 

S =Y g j dTdx l( d H dT ) 2 + (dcf/dx) 2 } - \N J drdxdytfir, x) ■ $(t, y)/\x - y\ a . (B2) 

We have set v = 1. g « 1/s is a coupling constant, not related to what we called g in other sections, g and A are 
scaled by N in order to have a smooth large- N limit. Inside a path integral, we may integrate over all fields 4>(t, x), 
without the constraint of Eq. I|B1|) provided that we introduce a Lagrange multiplier field, <t(t, x): 

S^S+— [ dTdx<j{$ 2 - 1). (B3) 
2ff J 

The action is now quadratic in unconstrained fields, so that we may, in principle, do the Gaussian integration over <j>. 
For this purpose it is convenient to write the long range term in uj-k space using: 

Here a is a short distance cut-off and this equation is valid for \k\a < 1. T is Euler's Gamma function. Note that 
the prefactor blows up, oc l/(a— 1) as a — ► 1. Note also that for 1 < a < 2, both T(2 — a) and sin[7r(2 — a)/2] are > 0, 
so that the second term in Eq. (52) is > 0. Asa-* 2, T(2 - a) -> 1/(2 - a) so that T(2 - a) sin[vr(2 - a)/2] -> tt/2. 
For 2 < a < 3, both T(2 — a) and sin[7r(2 — a)/2] < so that their product is positive, blowing up as a — ► 3. 

The first term in Eq. <|B4|I can be eliminated by shifting it by a (imaginary) constant. Thus we may write: 

A ' ' Lc ' ,k " w, k) ■ 0(-w, -k)[cj 2 + k 2 + CM^- 1 ] + ^ / drdxa($ 2 - 1). (B5) 



2g J (2tt)2^ ' ' ^ 1 1 ' 2g 

Here 

4 



C 



-r(2-a)sin[7r(2-a)/2] (B6) 
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and V is the space-time volume. Now we do the functional integral over 4>. This gives an effective action for the field 
a, which has an overall factor of N in front of it, and no further dependence on N. At large N the functional integral 
over a is dominated by a saddle point corresponding to a constant and purely imaginary value of a. Assuming a is 
constant, this effective action is: 

S ^' V = T {/ ISp + fc2 + c ^ + ^ - 7} < B7 > 

The saddle point is found by looking for a stationary point of S e ff- Setting ia — m 2 at the saddle point gives the 
self-consistent equation which determines m 2 : 



1 f dbjdk 
9 = J W 



g J (2tt) 2 lo 2 + k 2 + CX\k\ a - 1 +m 2 
rA dk 1 



-L 



4tt ^/k 2 + C\\k\ a ~ 1 +m 2 



(B8) 



Here we have introduced an ultra-violet cut-off, A of 0(a~ 1 ). For any a, A and g for which this equation has a solution 
with m 2 > 0, the system is in the disordered phase, with a finite gap, m. Note that, for a > 3, the integral diverges 
as k — > when m = 0. For small finite m it behaves as ln(A/m). Since this diverges as m — > 0, there will always be a 
solution for m, no matter how small is g. On the other hand, for 1 < a < 3, there will be a solution for A < A c only. 
At the critical value of A, m = 0, so A c is determined by: 



1 f A dk 1 



9 J -a 4tt ^ k 2 + CXc \k\^ 
As a — ► 1, C — > 4/(a — 1) and this becomes: 



1 /" A dfc 



5 i-A 4tt + 4A c /(a - 1) 



(B9) 



(BIO) 



This we see that A c oc (a — 1), as a — > 1. a = 1 is the critical value of a, for which A c — > 0. The behavior of 
A c (a) is qualitatively similar to what is obtained from SW theory, including the behavior near a = 1. Right at the 
critical point, the AF spin-correlation function is determined by the effective action with m = 0, and <f> treated as a 
non-interacting, free field, in the large N approximation. This implies a dispersion relation: 

uj = y/k 2 + CX\k\ a - 1 , (Bf f) 

and hence a dynamical exponent: 

z=(a-l)/2. (B12) 

This dispersion relation is the same as in SW theory, but there is no long range order at the critical point and hence 
no ambiguity in the value of z. The spin correlation function is given by: 

n f rlmrlh „i(uJT+kx) 

<'c*'w»>= p, l/isW' (B13) 

In particular, the equal-time correlation function is: 

<<j> a (0,x)<j> b {0,0)>=6 ab 4; I d \ ; ^ ■ (B14) 

V V >V K ' N J (4tt) ^k 2 + C\\k\^ V ' 

At large x, we may approximate this by dropping the k 2 term. It then follows from a rescaling that this decays as: 

< a (O, x)<p b (0, 0) >cx ^" a)/2 - (B15) 

The standard definition of the critical exponent rj then implies: 

z-l + ?7 = (3-a)/2, (B16) 
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giving: 



T] = 3 — a. 



(B17) 



Note that this is simply the behavior of the transverse correlation function in the ordered phase, according to SW 
theory. 

To calculate v, we need to calculate how m vanishes with A c — A as A — > A c , from below. It turns out that there are 
two different behaviors, depending on whether 1 < a < 5/3 or 5/3 < a < 3. A small change in A leads to a change in 
1/g which is linear in SX. Consider the effect of a small non-zero m on the gap equation. For 1 < a < 5/3, we may 
Taylor expand the gap equation: 



w 

9 J-. 



dk 
4ir 



(™ 2 /2)- 



(B18) 



(Note that this integral is finite at k — > for (3(a — l)/2 < 1 only, implying a < 5/3.) Setting A = A c , we see that: 

l/ 5 -l/.g c cx-m 2 . (B19) 
Thus to 2 oc — SX in this case. Now consider the equal time correlation function for a small m 2 . 



) y^k 2 + CX\k\ a - 1 + 



(B20) 



At large distances and small m 2 we should be able to drop the k 2 term and ignore the ultra-violet cut-off. A rescaling 
of the fc-integration variable then implies a correlation length: 



£ oc m- 2 ^- 1 ) oc (SX)- 1 /^ 

and hence an exponent: 

v = l/(a-l), (Ka<5/3). 
Now consider the case 5/3 < a < 3. Keeping a small non-zero m 2 , we find: 



dm 2 G) 2 / 



rffc 



1 



dm 2 \gj 2 7 47r[AC|fc| Q - 1 +m 2 ] 3 / 2 ' 



(B21) 
(B22) 

(B23) 



Since the integral is dominated by |fe| of 0(m 2 ^ a x \ we have taken the cut-off to infinity and dropped the k 2 term. 
By scaling, we see that: 



dm 2 \g 



oc m 



(5-3a)/(a-l) 



Integrating with respect to m 2 gives: 



for a constant, A. Thus we see that: 



This gives: 



Thus: 



I^l-Am^ 3 -^- 1 ), 
9 9c 



moc(-W)( a - 1 '/( 3 - a ). 
eocm-^-^ocC-^A)- 2 ^ 3 -"). 



3-a' 



(5/3 < a < 3). 



(B24) 

(B25) 

(B26) 
(B27) 

(B28) 
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